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^ ■ ABSTRACT 

o 

^ . The evolution of a gas shell, swept by the supernova remnant of a massive first 

. generation star, is studied with H2 and HD chemistry taken into account and with 

use of a semi-analytical approximation to the dynamics. When a first-generation 
^ . star, formed in a parent pregalactic cloud, explodes as a supernova with explosion 

^ . energy in the range of lO^^erg — 10^^ erg at redshifts of z = 10 — 50, H2 and HD 

molecules are formed in the swept gas shell at fractional abundances of ~ 10"'^ 
and ~ 10~^, respectively, and effectively cool the gas shell to temperatures of 
32K — 154 K. If the supernova remnant can sweep to gather the ambient gas of 
mass 6 x lO'^M© — 8 x 10^ Mq, the gas shell comes to be dominated by its self- 
gravity, and hence, is expected to fragment. The amount of swept gas necessary 
for fragmentation increases with the explosion energy and decreases with the 
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interstellar gas density (or redshift) of the host cloud, which provides a lower 
boundary to the mass of the host cloud in which star formation is triggered by the 
first-generation supernova. Also, the condition for fragmentation is very sensitive 
to the thermal state of interstellar gas. Our result shows that for a reasonable 
range of temperatures (200K ~ 1000 K) of interstellar gas, the formation of 
second-generation stars can be triggered by a single supernova or hypernova 
with explosion energy in the above range, in a primordial cloud of total (the dark 
and baryonic) mass as low as a few x 10^ Mq. For higher temperature in the 
interstellar gas, however, the condition for the fragmentation in the swept gas 
shell demands a larger supernova explosion energy. 

We also follow the subsequent contraction of the fragment pieces assuming 
their geometry (sphere and cylinder), and demonstrate that the Jeans masses 
in the fragments decrease to well below a solar mass by the time the fragments 
become optically thick to the H2 and HD lines. The fragments are then expected 
to break up into dense cores whose masses are comparable to the Jeans masses 
and collapse to form low mass stars that can survive to date. If the material 
in the gas shell is mixed well with the ejecta of the supernova, the shell and 
low-mass stars thus formed are likely to have metals of abundance [Fe/H] ~ —3 
on average. This metallicity is consistent with those of the extremely metal-poor 
stars found in the Galactic halo. Stars with such low metallicities of [Fe/H] < —5 
as HE0107-5240, recently discovered in the Galactic halo, are difficult to form by 
this mechanism, and must be produced in different situations. 

Subject headings: cosmology: theory — early universe — molecular processes — galaxies: 
formation — stars: formation 



1. Introduction 

Star formation in the early Universe is believed to have played a critical role in the 
formation and evolution of galaxies. Massive stars explode as supernovae that may contribute 
to the cosmic reionization and metal pollution of the Universe. Low-mass stars formed in the 
early Universe must survive to date, and are expected to carry precious information during, 
or even previous to, the early epoch of galaxy formation. Because of this importance, the 
process of star formation in the early Universe has been attracting wide interest (e.g., Ostriker 
& Gnedin 1996; Uehara et al. 1996). 

In the bottom-up scenarios of structure formation like those of the cold dark matter 
models, the first collapsed objects should have formed at redshifts of 2; ~ 10^ — 10 with 
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mass scales of 10^ Mq - 10^ Mq (e.g., Haiman et al. 1996; Tegmark et al. 1997). Stars 
ought to have been born in these first collapsed primordial clouds, totally lacking metals, 
before galaxies were formed. Such stars are referred to as first-generation stars. When 
the first collapsed objects are virialized, gas is once heated up to temperatures of T ^ 
100/1^/3 (M/lO^M0)2/3 (1 + z) K (Bromm, Coppi, & Larson 2002), and hence, needs to 
be cooled efficiently to collapse further into stars in these primordial clouds. Theoretical 
studies have suggested that in the first collapsed pregalactic objects, hydrogen molecules 
can form and cool the gas to temperatures of T ~ several x 10^ K, which leads to the 
formation of first-generation stars (e.g., Yoneyama 1972; Palla, Salpeter, & Stabler 1983). 
First-generation stars are expected to be massive, or low-mass deficient due to temperatures 
higher than those of the present interstellar matter (Bromm, Coppi, & Larson 1999; Bromm, 
Coppi, & Larson 2002; Nakamura & Umemura 1999; Abel, Bryan, & Norman 2000, but see 
Nakamura & Umemura 2001). 

On the other hand, the number of extremely metal-poor stars found in our Galaxy has 
recently been increased substantially (more than 100 for stars with [Fe/H] < —3) by large- 
scale surveys of the Galactic halo such as the HK survey (Beers et al. 1992 and see also 
Norris et al. 1999) and Hamburg/ESO survey (Christlieb et al. 2001). In particular, one 
giant star HE0107-5240 with a metallicity of [Fe/H] = —5.3 has been found very recently 
(Christlieb et al 2002). For such low metallicities as [Fe/H] ^ —4, metals contribute little to 
radiative coohng (Yoshii & Sabano 1980; Omukai 2000). The thermal property of such metal- 
poor gas has to be essentially the same as that of the primordial gas. The very existence 
of these low mass stars with such extremely low metallicities evidences the operation of the 
mechanism to form low-mass stars efficiently in the gas clouds of primordial abundances, 
completely devoid of metals. From existent studies on the thermal evolution of primordial 
gas clouds, the asserted lack in the low mass stars among the first-generation stars is mainly 
ascribed to the scantiness of relic electrons ([e/H] ~ 10^^ — 10~®), which limits the formation 
of H2 molecules necessary for cooling of gas in the primordial clouds (e.g. Galli & Palla 1998; 
Bromm, Coppi, & Larson 2002). The electron abundance can be larger, however, when gas 
in the primordial clouds is once heated above the temperatures of ^ 10^ K to be ionized and 
then undergoes coohng and recombination (Shapiro & Kang 1987). A larger fraction of free 
electrons will survive to yield production of more abundant H2 molecules primarily through 
H~ process of H -f e~ H^ + hu and H -|- H~ — > H2 + e^. Accordingly, the gas clouds 
can be cooled effectively first by H2 molecules and then by HD molecules to temperature 
sufficiently low for the formation of low mass stars. 

There are two possible ways to achieve such situations in primordial clouds (e.g. Shapiro 
& Kang 1987; Ferrara 1998; Nishi & Susa 1999; Uehara & Inutsuka 2000). For primordial 
clouds more massive than ~ IO^Mq, the virialized temperatures at the collapse can be 
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high enough to ionize the gas. For such massive clouds, however, we have to take into 
account the possible influence of pollution, prior to the collapse, with metals produced by 
the first-generation stars, since their collapse is delayed on average for more than ~ 10^ yr 
as compared with first collapsed clouds with masses of ~ IO^Mq (Tegmark et al. 1997). 
The other way may be shock heating by supernova explosions of the first-generation stars. 
The ambient gas in the clouds must be swept together, heated up and ionized by the shock 
waves induced by the supernova explosions. 

Tsujimoto, Shigeyama, & Yoshii (1999) have proposed a scenario of consecutive star 
formation, triggered by supcrnovac, from an abundance analysis of the extremely metal- 
poor stars found in the Galactic halo. It is to be properly established, however, that a 
gas shell swept by a supernova shock can actually fragment to form stars within the host 
clouds. Nishi & Susa (1999) discussed the condition of cloud disruption by supernovae with 
hydrogen chemistry taken into account. They argued that the scenario is possible only for 
massive first-collapsed objects with total masses larger than several xlO^M©. Since they 
did not solve the thermal evolution of supernova remnants, however, the electron density, 
and hence, the resultant H2 abundance may be underestimated, while the kinetic energy 
of supernova remnants is overestimated by applying the Sedov- Taylor solution to the later 
evolution. Ferrara (1998) discussed the notion that the gas, which is cooled to temperatures 
of ~ 300 K, is blown away from the cloud without fragmentation, in his study of the thermal 
evolution of a shocked gas, taking into account the collective supernovae of the Population 
III objects of total masses from 10^ Mq to 10^ Mq. 

Recently, Salvaterra, Ferrara, & Schneider (2003) have studied star formation induced 
by primordial supernova remnants (SNRs) in first collapsed objects, taking into account the 
H2 cooling. They suggested that a supernova shock can trigger formation of low-mass stars 
only when it is driven by an efficiently energetic supernova of explosion energy greater than 
lO^^erg (e.g., hypernova or pair- instability supernova). In their calculations, they assumed 
that the interstellar matter in which the SNR expands is ionized by UV radiation of the pro- 
genitor star and its temperature remains constant at lO'^K during expansion of the SNR. The 
resulting high interstellar pressure tends to stall the expansion of the SNR shell, preventing 
the shell from sweeping and gathering enough gas to become self-gravitating. However, the 
temperatures of interstellar gas are likely to decrease to lower values due to efficient radiative 
coohng by H2 molecules that are reformed in the interstellar gas on a timescale of less than 
10^ ~ 10^ yr after the supernova explosion. This reduction in the ambient pressure may 
significantly affect the expansion of the SNR shell and changes the condition to trigger star 
formation. Bromm, Yoshida, & Hernquist (2003) have also performed SPH simulations of 
such an SNR driven by a pair-instability supernova explosion due to a first-generation star. 
They showed that a high-energy supernova explosion (pair instability supernova) blows a 
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minihalo of mass 10^ Mq and gas with metallicity Z ^ 10~'^Zq is ejected to the intergalac- 
tic space. However, they did not address the issue of star formation triggered by the first 
supernova itself. 

In this previous work (Ferrara 1998; Salvaterra, Ferrara, & Schneider 2003; Bromm, 
Yoshida, & Hernquist 2003), only H2 molecules have been considered as the coolants after 
recombination. It should be pointed out, however, that under some conditions in the pri- 
mordial clouds, deuterated hydrogen molecules, HD, become a more efficient coolant than 
H2 (e.g., Puy et al. 1993; Galli & Palla 1998). For example, Uehara & Inutsuka (2000) 
have demonstrated that the HD cooling becomes dominant in the post-shock gas with shock 
velocities larger than ~ 300kms^^. Flower (2002) and Flower & Pineau des Forcts (2003) 
also investigated the effect of HD cooling on fragmentation of primordial clouds. Nakamura 
& Umemura (2002a) found that if H2 abundance exceeds a critical value of ~ 3 x 10~^, the 
thermal evolution of primordial clouds is controlled by HD cooling, which reduces the gas 
temperature to ~ 50 — 100 K (see also Nakamura & Umemura 2002b, in more detail). These 
studies indicate that HD cooling plays an important role in star formation in primordial gas 
clouds. 

In the present paper, we investigate the dynamical and thermal evolution of a gas 
shell swept up by an SNR in a first collapsed primordial cloud and discuss the possibility 
of low-mass star formation that can survive to date. We take into account the chemistry 
of H2 and HD molecules, and the evolution of an SNR is studied using a semi-analytic 
formula, which enables us to sweep the large parameter space with detailed treatment of the 
primordial chemistry. Uehara & Inutsuka (2000) have computed the evolution of a gas cloud 
for a particular initial condition, and demonstrated that once the gas has been ionized, the 
HD molecule plays a dominant role in cooling and enables low mass star formation. Our 
purpose is then to explore the conditions of explosion energy and the redshift of supernovae 
in primordial clouds that can lead to low mass star formation. First, we study the condition 
whether the gas shell, swept by an SNR, can break up into fragments as a result of cooling 
due to H2 and HD molecules. Further, we follow the evolution of the fragments into lower 
mass cores from which low-mass stars will be formed. In the computations, we also solve 
the thermal evolution of ambient gas simultaneously with ionization by UV radiation from a 
progenitor star, and discuss the effects of the thermal state of ambient gas on the evolution 
of the SNR. Since we are interested in low-mass star formation, triggered by a supernova 
explosion in a very first collapsed object, we work on the conditions appropriate to such 
situations. The approximation and method of computations are given in §2, and the results 
of our computations will be presented in §3. In §4, we derive the conditions under which 
supernova explosions of first-generation stars can trigger low-mass star formation in first 
collapsed clouds as functions of collapse redshift and energy of the explosion. We discuss the 



-6- 



relevance of our results in relation to the observed extremely metal-poor stars in §5. 



2. Model 

A massive first-generation star explodes as a supernova, sweeping up the ambient 
medium to form an expanding gas shell. If the gas shell becomes self-gravitating, it is 
likely to break up into fragments, which are expected to re-fragment into denser cores where 
the next-generation stars are formed. We explore the evolution of a supernova remnant 
(SNR) from the expansion of a shock wave in the ambient gas to the formation of dense 
cores in the fragments of the swept gas shell. The evolution of the SNR depends on two 
parameters: the explosion energy of the supernova, Eq, and the density of the interstellar 
gas, Pq, in the host cloud. For given parameters, wc solve the expansion and variations in 
the structure of the SNR by applying a semi-analytical approximation to the blast wave 
expanding in the uniform density. The chemical compositions in the gas shell swept by the 
supernova shock, which play a critical role in determining the thermal state of the SNR, 
are calculated by directly integrating the rate equations. To estimate the temperature and 
pressure in the ambient gas, we also solve the chemical and thermal evolution of the ambient 
gas. We describe our methods and approximations in the following subsections. 



2.1. Timescales 

We start by defining four typical timescales that characterize the evolution of an SNR, 
i.e., the expansion timescale Texp, the free-fall timescale ts, the dynamical timescale Tdyn, 
and the cooling timescale TcooI- 

The expansion timescale is defined as 

"^exp = -r , (1) 
rt 

where R and R denote the radius and the expansion speed of an SNR, respectively. The 
cooling timescale is given by 

" (7-1) A(T,n, Compositions) ' ^ ' 

where T and n are the temperature and the number density in the swept gas shell, and k 
and 7 are the Boltzmann constant and the ratio of specific heats, respectively. The symbol 
A (erg cm~^s~^) denotes the total energy loss rate per unit volume, which sums up all the 
cooling rates that are summarized in Appendix A. 
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The dynamical timescale represents the time in which the sound wave crosses the gas 
shell as 

_AR _ Rpo/3p 

where po is the density of the ambient gas; p, Cg, AT? and p denote the gas density, the sound 
speed in the shell, the shell width, and the mean molecular weight of the shell, respectively: 
ma is the atomic mass unit (p = pmau) . 

The free-fall timescale is written as 

where G is the gravitational constant and C is a structure parameter, given by C — 37r/32 
and C = S/tt for the spherical and cyhndrical collapses, respectively. Note that the free-fall 
timescale for the cylindrical collapse (ts^c) is larger than that for the spherical collapse (rff^g) 
by a factor of 2.94. See §2.2 for more detail. 



2.2. Expansion of SNR and Star Formation in the Fragments 

In our calculation, we treat the entire evolution from expansion of the supernova remnant 
shell, fragmentation of the swept gas shell and contraction of fragments in the expanding 
shell, through the formation of dense core consistently. In the following, wc divide the 
evolution into two phases, 'the ante-fragmentation phase' and 'the post- fragmentation phase^ 
for convenience's sake. The ante-fragmentation phase is defined as a period in which a 
supernova shock propagates and gathers interstellar gas until the swept gas shell is affected 
by gravitational instability and undergoes fragmentation. The post-fragmentation phase is 
defined as a period after fragmentation occurs until the fragments contract to be optically 
thick against the fine emissions by H2 and HD molecules. 



2.2.1. The ante- fragmentation phase 

Evolution of SNRs was studied extensively in the early 1970's using one-dimensional 
hydrodynamic codes (Cox 1972; Chevalier 1974). Summarizing their results, the evolution 
of SNRs is divided into three stages as (1) the free- expansion stage, (2) the Sedov-Taylor 
adiabatic stage, and (3) the pressure- driven expansion stage (see also Ostriker & McKee 
1988). As shown below, the expanding gas shell comes to be dominated by self-gravity and 
is expected to fragment during the pressure-driven expansion stage. 
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(1) the free-expansion stage (rexp <S Tcooi) 

In this stage, the SN ejecta expand freely. After the ambient gas swept by an SNR surpasses 
the SN ejecta in mass, the expansion is decelerated and the SNR enters the Sedov- Taylor 
adiabatic stage. With the ambient density po and the ejecta mass M^j, the transition radius 
is given by 

Ri = (3Mej/47rpo)'/'. (5) 

We put Mcj = 10 Mq, for simplicity. We have confirmed that the choice of the ejecta mass 
hardly affects the later evolution. 

From comparison with the Sedov- Taylor self-similar solution (see eq. [7] below), the 
transition time from the free-expansion to Sedov- Taylor stages is given as 



Ri 
1.15 



1/5- 



5/2 



1 



3Mej 



1.15 V47rpi 



1/3 



1/5- 



5/2 



(6) 



where Eq represents the SN explosion energy. 

(2) The Sedov-Taylor adiabatic stage (rexp < Tcooi) 

After the transition time ti, the SNR is dominated by the blast wave, and hence, the structure 
tends to a self-similar solution (Sedov 1946; Taylor 1950). The expansion of the shock front 
is well approximated to the Sedov-Taylor solution and described as 



it: = 1.15 (j-y't'/', 



(7) 



and the expansion speed of the SNR shock front is expressed by 



dt 



0.46( ^I'^'r^/s 



Po 



(8) 



This gives the post-shock pressure and shock temperature 



ps 



0.42 



Po 



Tr 



ps 



(7 + 1) \poJ 

8_/^\ / £o\ 



2/5 



25 k 



2/5 



PoJ 



as 



-6/5 



-6/5 



(9) 
(10) 



These equations are derived from the Rankine-Hugoniot relation and the Sedov-Taylor solu- 
tion, and indicate that the post-shock pressure and temperature decrease with the expansion 
of SNR (Ostriker Sz McKee 1988; see also Sakashita & Ikeuchi 1996). The post-shock gas is 
also cooled via radiative energy losses. Equation (2) enables us to estimate the cooling time 
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just behind the shock front. After the coohng time-scale becomes shorter than the expansion 
time-scale, i.e. TcooI < 7"exp, a cooled shell will form just inside the shock front and the SNR 
enters the pressure-driven expansion stage. We denote the time when TcooI decreases to be 
equal to Texp by ^2- 

(3) The pressure-driven expansion stage (tcooI < ''exp and Tdyn < Tg) 

The shell expansion is driven by the high pressure in the hot low-density cavity. The equation 
of motion of the shell is written as 

47rpo dR^R <? , ^ , ^ 

^ - 4rf (Pi, - Phc) , (11) 

where Pin and Phc mean the pressure in the inner cavity and the ambient gas pressure 
(Sakashita Sz Ikeuchi 1996), respectively. We assumed that the pressure inside the cavity 
decreases adiabatically as 

P,. = /^.(|)". (12) 

In this equation, R2 and P2 represent the radius of the shock front and the post-shock 
pressure, respectively, which are given from equations (7) and (9) at the beginning of the 
pressure-driven expansion stage {t = 12). The ambient gas pressure Phc is derived by calcu- 
lation of the chemical reaction and thermal evolution of ambient gas, as shown in §2.4. 

Using the above equation (11), we calculate the SNR expansion until fragmentation 
occurs in the shell (see below) or the expansion velocity of the SNR shell, V, becomes slower 
than the sound speed in the ambient gas of the host cloud, Cg^hc- As a result of the radiative 
cooling, the free-fall time scale tq grows shorter than the sound-crossing time (rdyn = AR/cg) 
over the shell width AR, and the self-gravity finally becomes dominant over the pressure 
force in the shell. Accordingly, we may well assume that the shell breaks into spherical 
fragments or cylindrical filaments when the following condition is satisfied: 

rs = rayn, or AP = c,(C/Gp)^/l (13) 

The shell width AR is related to the shell density p as Ai^ARR^p = (47r/3)P^po; where 
we used the mass conservation relation. Note that this gives AR = R/12 for the Sedov 
adiabatic stage if 7 = 5/3. As a corollary, the above condition = Tdyn is equivalent to the 
condition that the Jeans length in the shell is equal to the shell thickness. Such fragments 
produced by the break-up of the gas shell may correspond to the filamentary structures seen 
in the SNR, such as the Cygnus loop. Bromm, Yoshida, & Hernquist (2003) showed that the 
evolutions of SNRs of first-generation stars are well approximated by the analytical solutions 
used in the present paper. However, due to the restriction of numerical simulations (limited 
spatial resolutions) they could not explore the process after fragmentation. We solve the 
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contraction of the fragmented pieces with a semi- analytical method shown below and figure 
out the minimum mass of the second-generation stars. 



2.2.2. The post-fragmentation phase 

After the gas shell fragments, each fragment piece begins to contract in a free-fall 
timescale. In this phase, we consider two different geometries for the fragment: the spherical 
and cyhndrical configurations. For the spherical case, the contraction obeys the following 
equation: 

dv Grrir , 
dt " ^ ^ 

where v and represent the infall velocity and the mass contained within the radius r. For 
a uniform sphere, the timescale of contraction agrees with the free-fall timescale expressed 
by equation (4). Since the SNR shell fragments due to its self-gravity, on the other hand, the 
fragment is likely to have cylindrical geometry and the gas contracts in the radial direction, 
i.e., perpendicularly to the cylinder axis. Virial analysis of the cylindrical isothermal filament 
gives the equation of motion for the filament as follows: 

^ = -- (A - AJ (15) 

(Ostriker 1964; Uehara et al. 1996), where A denotes the line mass of the filament [= 
7r(Ai?/2)^p], and Ac is the critical value corresponding to the line mass of the filament in 
hydrostatic equilibrium and given by 

Ac = = 16) 

This equation indicates that the line mass density must be larger than Ac for the cylindrical 
fragment to contract. This condition is rewritten as 



in terms of the free-fall timescale for a cylinder with uniform density, tq^c — {^/'^GpY^'^ from 
equation (4). 

As shown in the next section, the main coolants in the fragments are H2 and HD 
line cooling. If the fragments are optically thick to these line emissions, the subsequent 
contraction proceeds nearly adiabatically, and along with the rise in the temperature, the 
right-hand side of equation (15) becomes no longer negative. Then, the radial contraction of 
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the fragment will be terminated, and the fragmentation is expected to occur again along the 
cylinder axis. Accordingly, we terminate following the contraction of the fragment when the 
fragments become optically thick to the line emissions of the most efficient coolant, either 
HD or H2, where the optical depth is evaluated under the escape probability method. In 
order to estimate the optical depths of H2 and HD lines, we use the Large Velocity Gradient 
(LVG) method (e.g., Goldreich & Kwan 1974). The optical depth for a transition J + 1 — > J 
is given by 



where h is the Plank constant, c is the speed of hght, Bj^j+i is the Einstein's B coefficient, 
and rij and gj arc the number density and statistical weight of the J-th level, respectively. 
The number density of the J-th level is evaluated using a two-level transition model with 
fitting formulas of the collisional de-excitation rates by Galli & Palla (1998) and Flower 
& Roueff (1999). The first six rotational transition levels arc taken into account. The 
velocity gradient, \dv/dr\, is approximated as avti^/ Rj = a{nGp)^^'^, where a, Vth, and Rj 
are non-dimensional numerical constant, the thermal velocity of HD, and the filament radius, 
respectively. The value of a depends on the velocity distribution of the collapsing filament. 
When the density reaches a critical density of HD cooling, the contraction of the filament 
becomes quasi-static and therefore, the velocity gradient becomes subsonic (Nakamura & 
Umemura 2002b). We thus set a to unity, for simplicity. 

It is thought that the Jeans mass, when the fragment becomes optically thick against 
the line emissions of the most efficient coolant (H2 or HD), gives the characteristic minimum 
mass for stars that will be formed in the fragment pieces, and hence, if it is much smaller 
than ~ 0.8 Mq, we may well expect that the low mass stars surviving today can be formed. 
For respective stages, we study the thermal and chemical histories of the SNR shell by a 
one-zone approximation, which is explained below. 



The temperature at the transition epoch, ti, from the free-expansion to the Sedov- Taylor 
adiabatic stage, is determined from the condition that the transition is continuous as 




(18) 



2.3. Temperature and Density in the Shell 



2.3.1. The ante- fragmentation phase 



(19) 
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with 7 = 5/3. We calculate the variation in the temperature of the shell under the one- zone 
approximation in the following way: 

dT fdT\ fdT\ fdT\ 

where {dT/dt)e^p, {dT/dt)y.a,d, and {dT/dt)comp represent the term due to the expansion cool- 
ing, the radiative cooling, and the compressional heating due to contraction of the shell, 
respectively. In the Sedov-Taylor adiabatic stage, the postshock temperature is given by the 
similarity solution, where T oc t~^/^. Thus, the expansion cooling term is denoted as 

y„r-'"'- '''' 

This term is only effective in the Sedov-Taylor adiabatic stage, and in the pressure-driven 
expansion stage, tends to be negligible as compared with the radiative cooling and the 
compressional heating term. Thus, we include the expansion cooling only in the Sedov- 
Taylor stage and ignore it after the pressure-driven stage. 

The radiative coohng term is denoted as 

( c(^\ ^ ^ ^ y (7 - 1) A(T, n. Composition) 

In this cooling, we include the inverse Compton cooling (Ikeuchi & Ostrikcr 1986) and 
the radiative cooling by the atoms and ions of H and He (Cen 1992) and by H2 and HD 
molecules (Galli & Palla 1998; Flower et al. 2000). We take into account the effect of the 
CMB radiation on the radiative cooling by modifying the total coohng rate as 

A = A(r)-A(TcMB), (23) 

where the cosmic background temperature Tcmb is calculated as Tcmb = 2.73(1-1-2;) K. 
Owing to this formula, the gas cannot cool below the CMB temperature at that time. 

The compressional heating has a form 

dT\ . .T dn 

where n represents the number density of molecules, atoms, ions and electrons. This term 
becomes effective after the shell has cooled and the density increases. 

The density in the shell is given by the jump condition across the strong shock front 
during the Sedov-Taylor stage as 

7 1 

P - ^Po, (25) 



where po represents the ambient gas density. This equation indicates that the density in 
the shell is equal to 4po for 7 = 5/3. After the condition TcooI < T"exp is realized in the 
pressure- driven stage, time evolution of physical quantities in the shell differs according to 
whether the cooling time is shorter than the sound crossing time of the shell (tcqoI < Tdyn) 
or the other way around (xdyn < Tcooi)- In the former case, the shell cools without changing 
the width of the shell. Thus, we assume 

^ = 0, for Tcool < Tdyn- (26) 

Accordingly, the pressure in the shell varies in proportion to the temperature given by 
equation (20) in this regime. 

If the sound-crossing time is shorter than the cooling time (rdyn < Tcooi), on the contrary, 

the shell is subject to compression, and in the decelerating frame with the gas shell, the 
structure tends to be in hydrostatic equilibrium with the boundary pressures given by the 
inner cavity pressure (eq. [12]) and the post-shock pressure. Pout- Thus, we assume that the 
pressure in the shell is brought near to the average of Pin and Pout in a dynamical timescale 
as follows: 

dP /^Pin + Pout ^\ 1 



, o ~ ^ ' fo'^ '^dyn < Tcool, (27) 

at \ 2 J Tdyn 

where the post-shock pressure is assumed as the sum of the ram pressure and the pressure 
in the ambient gas Phc: 

Pout = ^-rPoP' + Phc- (28) 
7 + 1 

The density variation at this stage is determined by equations (20) and (27) with the equation 
of state. 



2.3.2. The post-fragmentation phase 

During the pressure-driven expansion stage, the dynamical timescale lengthens as the 
shell expands, eventually reaching the free-fall timescale. When the condition Tg < Tdyn 
is satisfied, the shell begins to fragment owing to self-gravity, and the evolution enters the 
final stage of the fragmentation and star formation stage. As mentioned above, the subse- 
quent evolution of the fragments seems to differ according to their geometry. We consider 
two different types of fragment geometry, i.e., the spherical and the cylindrical symmetric 
configurations. For spherical fragments, the time variation of the density is approximated as 
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with the free-fall timescale given by equation (4) . 

On the other hand, if the fragments form long cyhndrical shapes, the widths of the 
filaments shrink with time according to equation (15). Since equation (15) gives the contrac- 
tion time scale of ~ r/[2G(A — Xc)]^^"^ — l/[27rGp(l — Ac/A)]-^/^, we assume that the density 
changes as 

" ^-^V". (A>A.) (30) 



with the line mass density A determined at the fragmentation (Ostriker 1964). 

The cylindrical collapse proceeds much slower than that of the spherical one (Uehara et 
al. 1996). This comes from the factor 1 — Ac/A in equations (15) and (30). Along with the 
cooling of the gas, Ac decreases. If the critical line mass Ac becomes smaller than the line 
mass, i.e., A > Ac, the filament contracts. At the same time, the contraction causes Ac to 
increase because of the compressional heating. If Ac increases to approach A, the contraction 
slows down, and finally halts when Ac reaches A. Accordingly the contraction of filaments 
proceeds in the cooling timescale, keeping the temperature slightly lower than determined 
from the condition Ac = A, and hence, much more slowly than in the free-fall timescale. This 
makes a clear contrast to the fast collapse accomphshed for the spherical system [eq.(29)]. 

In our calculation, the initial line mass A is determined at the fragmentation epoch 
(rg = Tdyn) under the assumption that fragmentation occurs with a wavelength equal to the 
shell width. This gives a fine mass of 

where we used the shell thickness of AR — poR/3p. As the shell expands, the fragments 
accumulate additional gas that is swept by the shell. We assume that the line mass increases 
in proportion to the shell mass as 

dX _ 



where Ms„ and Mg^ represent the shell mass and the rate of increases in the shell mass, 
respectively. Along with the increase in A, the contraction becomes shorter than the expan- 
sion, and thereafter, A stays constant, since the contraction is further accelerated with the 
decrease in the free-fall and coohng timescales with the increases in the density. We also 
assume that the mass of the spherical fragments increase in proportion to the shell mass. 
However, since the spherical collapse is much faster than the cylindrical one, this does not 
play an important role. 



- 15 - 



2.4. The evolution of the ambient gas 

The thermal state of the ambient gas may play a part in the evolution of the expanding 
gas shell in the pressure-driven expansion phase (§2.2.1 (3)) and the resultant fragmentation 
of the swept shell. In the course of formation of first-generation stars, gas in the primordial 
clouds suffers cooling by H2 molecules formed with relic electrons as agents, and the gas 
temperature can decrease as low as ~ 200 K (e.g., Bromm, Coppi, & Larson 2002). Once 
a massive star is born, however, gas is ionized by the radiation of the massive progenitor 
stars and heated to T ~ lO^K. After the massive star explodes, radiation coming from the 
interior of the SNR may heat the ambient gas or destroy the hydrogen molecules. 

As for the UV radiation from a progenitor star, the extent of the ionized region depends 
on the ratio of the total number of ionization photons per unit time to the recombination 
rate of the hydrogen. If we assume an 05-type progenitor star and the ambient gas density 
of no = 1 cm~^, this star ionizes the ambient gas inside the Strongren radius of ~ 100 
pc (Panagia 1973). After the massive star explodes, the ambient ionized gas begins to 
recombine. Although the ambient gas continues to cool by radiative cooling, the gas near 
the shock front (preshock gas) may be heated by the ionization radiation from the SNR's 
hot interior (Shelton 1999, Slavin et al. 2000). 

The preshock gas ionized by the radiation emitted from the postshock hot gas is cal- 
culated by ShuU & McKee (1979) for a metal- rich gas cloud and by ShuU & Silk (1979) 
for a metal-poor gas cloud. The controlling parameter for this problem is the ratio of the 
emergent photon flux $(cm~^s^^) from the postshock gas to the hydrogen atom number flux 
noV{cm~^s''^) flowing into the SNR. ShuU & Silk (1979) showed that the emitted photons 
scarcely ionize the ambient gas for no — 1 cm~^ if the shock velocity is as low as F ^ 60 
kms~^, because the above photon-to-gas ratio ^/noV ^ 1 is too small at that time. Thus, 
the heating of the preshock gas by the SNR shock is negligible in the late evolutionary phase 
{V <^ 60km s^^). In the early evolutionary phase, the emission from the SNR shock is also 
negligible, because the shock is strong and the ambient gas pressure does not play a role, as 
discussed in Shapiro & Kang (1987). 

In our calculation, as long as the shock speed is faster than V > 60 kms~^, the ionization 
level of the ambient gas is always much higher than that expected from the ionization by 
the SNR (Shull & McKee 1979) and its temperature is maintained at The ~ lO'^K, because 
the ambient gas photoionized by the progenitor star has not recombined/cooled sufficiently 
by this epoch. Thus, we expect that SNR evolution is not affected, whether the ionizing 
photons emitted from the SNR are considered or not. Therefore, we neglect ionization by 
the ionization photons from the SNR. 
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The SNR radiation may prevent the ambient gas from cooling through the dissociation 
of the H2 molecules, which is the effective coolant at lower temperature (T < lO^K). The 
dissociation timcscalc of the molecular hydrogen, t^is, is given by ^dis — 2.8 x 10~^^F£^ yr, 
(Omukai & Nishi 1999). Here Flw (erg s~^ cm~^ Hz~^) is the average radiation flux in the 
Lyman - Werner (LW) bands. The flux Flw is estimated as a function of the shock velocity 
of SNR (ShuU & McKee 1979; ShuU & Silk 1979). The hydrogen molecules begin to form 
after the shock velocity decelerates below V — 60 km in our calculation. Prom the flux 
Flw — 1.3 X 10"^^ erg s~^ cm~^ Hz~^ for no = 1 cm~^ at the stage of F = 60 km s~^ 
(ShuU & Silk 1979), the dissociation timescale is given as t^is — 2.2 x 10^ yr, which is much 
longer than the formation timescale of the molecule hydrogen, nB.2/idnB.2/dt). After that, 
the dissociation timescale continues to increase, because the flux Flw decreases with time, 
and keeps longer than the formation timcscalc. Therefore, radiation from an SNR barely 
decreases the molecular fraction of ambient gas. 

As a result, we can safely ignore the effect of radiation from hot gas in the SNR and 
include only the effects of photoionization and heating due to the progenitor star. How- 
ever, the above effect may delay cooling of the ambient gas slightly. In §4, we discuss how 
the expansion of the shell is affected if we assume a higher-temperature ambient gas than 
calculated in the next section §3, and show the condition for low mass star formation. 

In order to obtain the proper boundary conditions of SNR shell evolution, we need 
to know the evolution of the ambient pressure. For this reason, we solve the chemical 
reactions and the thermal evolution of the ambient gas, simultaneously with the dynamical 
and thermal evolution of the SNR shell. The initial chemical composition of the ambient 
gas is derived under the assumption that the ambient gas has been heated and is kept at a 
temperature of The = 10^ K by the progenitor star of the first generation which lives for 10^ 
yr. After the progenitor star explodes, we assume that the heating source disappears and 
follow the cooling of the ambient gas by solving the same rate equations and the equation of 
energy conservation, as described above, under the constant density of po- Prom the ambient 
temperature. The, we derive the interstellar pressure Phc by 

-Phc = — ^- — The, (33) 
/ihc rria 

where /^hc is the mean molecular weight of the ambient gas. This pressure is applied to 
the evolution of the SNR gas shell in cq. (11). The SNR shell can continue expanding and 
sweeping the ambient gas, only when the expansion speed is faster than the sound speed as 

V > c,,hc- (34) 

Otherwise {V < Cs,hc), the gas shell dissolves and merges into the ambient gas. Accordingly, 
the above condition (eq.[34]) has to be satisfied for the gas shell to fragment. 
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2.5. Chemical Composition 

We have modeled the thermal and dynamical evolution of the SNR shell and ambient 
medium under the one-zone approximation, as stated above. The cooling rates in the gas 
shell and ambient medium are determined by their chemical composition. To estimate the 
abundances of chemical species in these gases, we solve the time-dependent chemical reaction 
equations numerically. In this paper, we consider the chemical reactions of the following 12 
species, H, H+, H~, He, He"*", He"*""*", H2, D, D"*", HD, HD+, and e~. We adopt the primordial 
composition obtained by Galli & Palla (1998) as the initial condition of our calculations. 
The reaction rates we include in our calculations are summarized in Appendix B. 

3. Results 

We begin our computation from the transition epoch, ti, between the free-expansion 
and the Sedov- Taylor adiabatic stages given by equation (6) for a given set of values of the 
explosion energy, Sq, and the density, po, of the host cloud. Then, we solve the equations 
of structural changes in the expanding SNR and the rate equations for the changes in the 
chemical abundances, formulated in the preceding section. In the present work, we adopt 
an explosion energy between sq = lO^^erg and lO^^erg, ranging from a normal supernova 
to a hypernova (Nomoto et al. 1999). The density of interstellar gas in the host cloud is 
taken to be in the range between po — 3.51 x 10~^^gcm^^ and 3.49 x 10~^^gcm~^ taking 
into account the formation epoch of the host clouds. Since the densities in the virialized 
host clouds are approximated to 200 times the average baryon density in the Universe at 
that epoch (White et al. 1993), the above density range corresponds to the redshifts of 
formation epoch 5; = 10 — 50 (po oc [1 + z]^) for the Einstein-de-Sitter Universe model with 
a Hubble constant of Hq — 70 km s~^ Mpc~^ and the baryon fraction of Ob — 0.06. 

The model parameters adopted are hsted in Table 1 with the characteristic physical 
quantities of gas shells and ambient gas when the condition of fragmentation, tq — Tdyn, 
is fulfilled. In the following, we discuss the results, dividing the evolution into two phases, 
defined in the preceding section; the ante-fragmentation phase (rg > r^yn) until the gas shell 
undergoes fragmentation and the post-fragmentation phase (rg < r^yn) during which the 
fragments collapse gravitationally, eventually forming second-generation stars. 
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3.1. The Ante- Fragmentation Phase 

Figure 1 shows the evolution of the shell radius and the expansion velocity of the SNR 
against the elapsed time from the SN explosion for various models. For each model, the 
expansion of the SNR is well expressed by two power-laws, in which the break occurs at the 
transition time t2 from the Sedov-Taylor (R oc i^/^; eq. [7]) to the pressure-driven stages 
{R oc i^/^; eq. [11]). The expansion velocity, and hence, the radius of the SNR, are larger 
for smaller collapse redshifts and for larger explosion energy. For example, the initial radius 
at the transition epoch ti = 2.1 x lO^yr for z = 10 is 4.7 times as large as that at ti — 
2.6 X lO^yr for z = 50. We note that the expansion velocity and radius depend weakly on 
ambient density, as inferred from the Sedov-Taylor similarity solution [R oc (eq/ PqY^^], but 
the ambient density depends strongly on the collapse redshift [p oc (1 + z)'^]. Therefore, the 
expansion velocity and radius change greatly with z. More quantitative results of all models 
are summarized in Table 1. 

The evolutionary variations of thermal state in the swept gas shell and ambient gas are 
shown in Figures 2, 3, and 4 as a function of elapsed time. Figure 2 shows the temperature 
variations of gas shell (top and middle panels) and ambient gas (bottom panel) for models 
with different redshifts (top and bottom panels) and for different explosion energies (middle 
panel) . The temperature of the SNR shell (top and middle panels) decreases gradually with 
expansion during the Sedov-Taylor adiabatic stage. The initial temperature is determined 
from the energy conservation at transition time ti in equation (19). The temperature steeply 
descends from T ~ 10^ K to lO'' K owing to efficient atomic cooling by He and H (see Fig. 7 
left panel). As the temperature decreases below 10^ K, the temperature drop slows again, 
and yet, continues to decrease below 100 K owing to cooling by H2 molecules and then by 
HD molecules. There are two small dents discernible around T ~ 10^ K and 150 K, which 
are caused by the coohng due to H2 and HD molecules, respectively. For higher redshifts, 
the temperature decreases to the CMB temperature. For £0 = 10^^ erg models, this occurs 
for z ^ 30, and then, the SNRs evolve isothermally, as seen in the top panel of Fig. 2 (at 
t ^ 10^ yr). For the larger explosion energy, however, the gas shell sweeps the ambient gas 
at a higher rate because of larger expansion velocity, and hence, can satisfy the condition 
for fragmentation at higher temperatures before it cools to the CMB temperature even for 
the largest redshift ^ = 50 as seen from the Sq> 1 x 10^^ erg models. 

Since the radiative cooling is in proportion to the density squared, the cooling is faster 
for models with larger ambient densities (or larger redshifts). Accordingly, the time interval 
from the supernova explosion to the shell fragmentation epoch is shorter for larger ambient 
gas density or for larger redshifts of the host cloud collapse, which decreases by a factor of 
10 — 30 from several x 10^ yr to ~ 10^ yr between z = 10 and 50. If the ambient gas density 
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is the same, the model with smaller explosion energy onsets the rapid atomic cooling earlier. 
After the gas cools to T ^ lO^K, the evolution converges because of the strong temperature 
dependence on atomic cooling. Accordingly, the time interval from supernova explosion to 
fragmentation is elongated for a smaller explosion energy, since it takes longer to gather the 
necessary mass because of smaller expansion velocity. 

In the bottom panel of Figure 2, we show the evolutions of the temperature, T^c, in 
the ambient gas for several different values of collapse redshift z. The ambient gas evolves 
almost isothermally for t ~ 4 x 10"'' ~ 10*^ yrs, and then, the temperature decreases below 
The = 1000 K for 10® ~ 10^ yr. For the model with z = 20 and Eq = lO'^'^erg, the temperature 
of the ambient gas decreases to The = 223 K in 14.7 Myr by the time the temperature of the 
SNR shell decreases to T = 53 K and satisfies the conditions of fragmentation. The cooling 
of the ambient gas is also faster for models with higher redshifts because of higher density of 
ambient gas. The ambient gas can finally cool to the temperatures The = 193 — 656 K when 
the condition of the fragmentation is satisfied, as listed in the 11-th column of Table 1 with 
the sound speed, Cs,hc, in the 12-th column. 

Figure 3 illustrates the time variations in the pressure in the gas shell, P, along with 
those of the boundary pressures, i.e., the cavity pressure inside the gas shell. Pin (eq. [12]) 
and the post-shock pressure outside the gas shell Pout (eq. [28]) for the model with z = 20 
and ^0 = 10^^ erg. As the atomic cooling becomes very effective in the pressure-driven 
stage, the pressure in the shell starts to decrease nearly in proportion to the temperature in 
the gas shell since TcqoI ^ Tdyn and hence, the dynamical readjustment of structure cannot 
catch up with the temperature drop caused by the radiative cooling. The radiative cooling 
rate declines very steeply below 10'^ K. As the cooling timescale grows longer and eventually 
exceeds the dynamical timescale (rdyn ^ Tcooi), the gas shell starts to readjust its structure to 
the surrounding pressures, and the pressure approaches the average of Pjn and Pout according 
to equation (27). In particular, for t ^ 10® yr, when the time grows longer than Tdyn, the 
gas shell restores the adjustment and the pressure P tends to be controlled by the boundary 
pressures and decreases with them. 

In this figure, we also plot the ambient pressure Phc- The ambient pressure remains much 
lower than the inner pressure, Pjn, for most of the time, and it is only after t > 1.3 x lO^yr 
that the ambient pressure grows higher than the inner pressure. Since their difference attains 
only 9 % at the time of fragmentation and is much smaller than the momentum flux of the gas 
shell, the effect of decelerating the expansion of the SNR shells is rather small, and hence, will 
not affect the evolution of the SNR shell so much. In this model, the SNR gas shell keeps 
expanding at a velocity exceeding the sound speed, and hence, undergoes fragmentation 
without dissolving and merging into the ambient gas. For the models of larger explosion 
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energy and higher ambient density, the condition that V > Cs,hc holds by a greater margin 
because of greater expansion velocity at the time of fragmentation, as seen from Table 1. 

The changes in the number density in the gas shell are plotted in Figure 4. The gas 
density in the shell is kept nearly constant at 4po until Tdyn becomes shorter than Tcooi at 
i ;^ 3 X 10^ yr and for T ^ 10^ K. Then, the density starts to increase as the gas shell 
undergoes compression by the pressures at the inner and outer boundaries. Finally, when 
the gas shell restores the hydrostatic equilibrium, the density begins to decrease along with 
the boundary pressures that confine the gas shell, as seen from the model with the smallest 
explosion energy. For the greater explosion energy, however, this final stage does not occur; 
because of greater expansion velocity, the SNR gathers more gas to make the shell thicker 
and the dynamical timescale longer in proportion, and hence, the fragmentation condition 
is satisfied before the hydrostatic equilibrium is restored. 

We exemplify the behaviors of the four timescales, mentioned above, in Figure 5, which 
shows the evolutionary changes in the expansion time, the cooling time, the dynamical time, 
and the free-fall time, mentioned above for three models of {z, Eq) = (20, 10^^ erg; left panel), 
(50, 10^^ erg; top- right panel) and (20, 10^^ erg; bottom-right panel). As the temperature 
descends to T ~ 10^ K, TcooI starts to decline rapidly owing to atomic cooling, and eventually 
becomes shorter than the expansion timescale Tcxp which increases with time. This causes 
a transition from the Sedov- Taylor adiabatic stage to the pressure-driven stage. When the 
temperature falls further to T ~ 2 x 10^ K, TcqoI reverses the direction of the change and 
begins to increase. As the recombination of hydrogen atoms proceeds with a further drop 
in the temperature, TcooI grows large very rapidly, increasing by a factor of ~ 100 while 
the temperature decreases to T ~ 5 x 10^ K, and forms a sharp minimum. Thereafter, 
the cooling timescale continues to increase constantly, slowing the increase rate slightly by 
the enhancement of coohng by H2 and HD molecules (two small dents are discernible at 
i ~ 5 X lO^yr and 5 x 10^ yr in the left panel). On the other hand, the dynamical timescale 
grows large owing both to the decrease in the temperature and to the increase in the shell 
thickness. When the atomic cooling begins to be effective, therefore, the cooling timescale 
becomes smaller than the dynamical timescale. As the temperature decreases below 10^ 
K, they reverse their relation and the cooling timescale grows larger than the dynamical 
timescale again. When Tdyn ^ Tcoob the gas shell undergoes compression, as seen above. 
Then the density in the shell increases, which causes a decrease in the free-fall timescale, 
as seen from i 5 x 10^ yr, while the dynamical timescale continues to increase because of 
the decrease in the shell temperature, and also because of the increase in the shell thickness. 
Finally, ts becomes as short as Tdyn {t — 10^ yr for the left panel), and the gas shell is 
expected to break up into fragments, each of which undergoes gravitational contraction. 



- 21 - 



For higher ambient density (or higher z] top- right panel), the coohng, dynamical, and 
free-fall timescales are shorter, while the expansion timescale is slightly dependent on the 
density in the gas shell. These accelerate the evolution. Accordingly, the fragmentation epoch 
becomes earlier, because the gas in the shell cools more rapidly and the free-fall timescale is 
smaller for higher density. For the model with greater explosion energy (bottom-right panel), 
the temperature remains higher, which delays transition to the pressure-driven stage. In the 
later stage, the difference in the temperature tends to be eliminated as a result of the atomic 
cooling. On the other hand, faster expansion makes the shell thicker and the dynamical 
timescale longer. This delays the arrival of the epoch of TcooI > 7"dyn and the contraction of 
the gas shell, as seen in Figure 4, but brings forward the epoch of fragmentation shghtly. 
Accordingly, we conclude that the time necessary for fragmentation depends significantly on 
the ambient gas density (po or the formation redshift z), but weakly on the explosion energy 

We present typical patterns of the evolution of fractional abundances and of the vari- 
ations in the cooling rates against the shell temperature, T, in the left panel of Figures 6 
and 7, respectively, for the {z, Sq) = (20, 10^^ erg) model. The ambient gas is ionized by the 
SNR shock at first, and then, recombines around T ~ 5 x 10^ K. Thereafter, H2 and HD 
molecules are formed through the reactions as 

e- + H ^ H- + hi/, (35) 

H- + H^H2 + e-, (36) 

and 

D+ + H2 ^ HD + H+. (37) 

Since H~ ions react as a catalyst to make hydrogen molecules, both H2 and HD abundances 
increase while H^ abundance increases, as in the left panel of Figure 6. The HD abundance 
continues to increase even after the H~ abundance decreases and the H2 abundance saturates. 
This is due to the charge transfer reaction in equation (37), through which HD molecules 
are formed from much more abundant H2 molecules. 

The left panel of Figure 7 plots the contributions from the cooling rates A (T, n, Composition) 
for the same model as in Figure 6. This elucidates the fact that H2 and HD molecules are 
effective coolants below the temperature ~ 2 x 10^ K and ~ 150 K, respectively. The cooling 
rate is smaller by a factor of 10^ and 10^ than the atomic cooling rate at its peak. If it 
were not for H2 and HD molecules, however, gas could not be cooled below ^ 10^ K within 
the Hubble time. In the model presented in Figure 7, the fractions of H2 and HD molecules 
amount to n(H2)/no,H = 1-58 x 10~^ and n(HD)/no,H = 8.32 x 10~^ at the fragmentation 
epoch, which are much larger than their primordial values of 10~^ and 10~^, respectively 
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(Galli &: Palla 1998). These abundances agree well with the result of Uehara & Inutsuka 
(2000). When the temperature becomes as low as 145 K, the HD cooling exceeds the H2 cool- 
ing and promotes a further drop in gas temperature. The fragmentation condition, = tdyn, 
is satisfied only when the temperature descends to T = 53 K at the age of 14.7 Myr after 
the SN explosion. By this stage, the shell expands to a radius of 90 pc and the gas shell has 
gathered the ambient gas of mass of 1.11 x 10^ Mq for this particular model. 



3.2. The Post-Fragmentation Phase 

After the condition of = Tdyn is fulfilled, we assume fragmentation to occur. Then 
the evolution proceeds to the final stage of fragmentation and star formation. We follow the 
gravitational contraction of the fragment pieces assuming two different types of fragment 
geometry: sphere and cylinder. 

The temperature evolution during the post-fragmentation phase of the{z, Eq) = (20, 10^^ 
erg) model is shown in Figure 8 for both the spherical and cylindrical collapses. For the 
cylindrical case, the gas temperature Tc decreases nearly to, but slightly higher than, the 
CMB temperature. Fragmentation occurs at the age of 14.7 Myr after the SN explosion. 
At this epoch, ~ 53K and ric ^ 7cm~^, which corresponds to a steep break near the 
lower-left corner in Figure 8. Before the fragmentation, the evolutionary path in Figure 8 
is vertically downward, that is, gas cools without strong compression. After this epoch, the 
fragment changes its evolution as it evolves nearly isothermally. The temperature of the 
fragment decreases further to 48K, which decreases the critical line-mass Ac oc Tc by 19% in 
6.7 Myr from the fragmentation epoch. The shell expands from ~ 90 pc at t = 14.7 Myr 
to i? ~ 103 pc at t = 28.7 Myr, which increases the line-mass A approximately 50%. These 
two factors work cooperatively to promote further collapse of the fragment. Thereafter, the 
fragment evolves under the condition that Ac ~ A. 

The gas temperature remains nearly constant for the number density n <^ 10^ cm^^, 
and then, increases to ~ 100 K owing to the increase in the mean molecular weight, 
as also seen from Uehara & Inutsuka (2000). Variations in the fractional abundances are 
shown in the right panel of Figure 6 against the number density in the fragment pieces 
for the cylindrical contraction case. As the fragment contracts, the HD fraction gradually 
increases with density through the reaction in equation (37). The H2 fraction changes little 
until the density becomes sufficiently high for the three-body reactions (n ^ lO^cm"^). A 
break near n ^ lO^cm"'^ in Figure 8 corresponds to this critical density of H2. Beyond that 
density, the H2 fraction starts to increase rapidly with density. The variations in the cooling 
rates are also plotted in the right panels of Figure 7 against the number density for both 
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the spherical and cylindrical cases. H2 and HD molecules are only the coolants effective 
after the fragmentation. It is clearly shown that HD molecules play a more important 
role in the lower temperature. HD molecules continue to dominate the cooling throughout 
the post-fragmentation phase, which allows the contraction of the fragments at much lower 
temperatures than found in the previous computations involving only H2 molecules (Puy et 
al. 1993; Uehara et al. 1996; Nakamura & Umemura 1999; Flower & Pineau des Forets 
2003). 

For the spherical case, on the other hand, the temperature keeps rising with density 
owing to dynamical compression (Fig. 8). When the temperature becomes higher than Tg ~ 
150K, HD molecules are dissociated to decrease the fractional abundance by a factor of 
~ 1/10, and hence, the coohng rate due to HD molecules is overtaken by that due to H2 
molecules for n > 10^ cm~^. Thereafter, the H2 fraction increases through the three-body 
reactions, which enhances the contribution of H2 molecules to the cooling. Although the HD 
fraction also increases again through the charge exchange reaction with the H2 molecules, the 
latter continues to dominate the cooling because of high temperatures. The temperature in 
the fragment pieces is, however, kept much lower than obtained in the previous computations 
without the HD molecules because of the lower temperatures at the fragmentation epoch 
(Ferrara 1998). 

We also show the time variations in the Jeans masses in fragment pieces during the post- 
fragmentation phase in Figure 8. For the cylindrical case, the Jeans mass Mjc decreases 
to 0.16 M0 at the stage when fragment pieces become optically thick at the density n ~ 
10^°cm~^. After this stage, the cooling becomes ineffective and the gas cannot contract 
in the mass-scale smaller than the Jeans mass, while the gas element more massive than 
the Jeans mass may collapse in a free-fall timescale. Because an adiabatic core is formed, 
further fragmentation can give rise to stars with the Jeans mass at the core formation epoch, 
which is well below one solar mass and corresponds to a less massive star which can survive 
to date. We note that our Jeans mass is larger than that of Uehara & Inutsuka (2000), 
Mjc ~ 0.04 Mq , because they seem to underestimate the optical depths of HD lines (see 
e.g., Nakamura & Umemura 2002). Therefore, our fragments are likely to evolve into low- 
mass stars, rather than primordial brown dwarfs. For the spherical case, on the contrary, 
the fragment may not stop contracting even after it becomes optically thick (e.g. Omukai 
& Nishi 1998). However, since the gas shell is thin (the ratio of the thickness to the radius 
~ 0.1), the fragment pieces are expected to form a cylindrical filament whose axis is parallel 
to the shell, and hence, can give rise to stars with sub-solar masses. 
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4. Fragmentation Conditions in the Primordial Clouds 

In the preceding section, we have solved the evolution of gas shells swept by the super- 
nova remnants of first-generation stars, assuming that the host clouds have sufficiently large 
masses. In order for star formation to be actually triggered in host clouds, the following two 
conditions have to be satisfied at the fragmentation epoch: i.e., (1) the shell expansion ve- 
locity is larger than the sound speed of ambient gas and (2) the mass of swept gas is smaller 
than the total baryon mass of the host cloud. In our calculations, we have shown that the 
SNRs can give rise to low-mass star formation while the expansion velocity is larger than 
the ambient sound speed, and hence, without dissolving into the interstellar gas, except for 
the two cases of weakest explosion energy of eg = 1 ^-^d 3 x 10^^ erg at the lowest collapse 
redshift of z = 10 (see V and Cs,hc in the 7th and 12th columns in Table 1). Then, the 
fragmentation occurs in the host cloud, if the baryon mass of the host cloud is larger than 
the mass, Mg^, of the gas swept by a SNR at the fragmentation epoch. If the expansion 
velocity is slowed to the ambient sound speed, the gas shell dissolves into the ambient gas. If 
a sufficient mass of gas is not available, SNRs expand beyond the edges of host clouds, either 
falling back and being mixed into the interstellar clouds or being dispersed out of clouds and 
spread into the intercloud space, depending on whether the expansion velocity is lower or 
greater than the escape velocity of the host clouds. We discuss below first the dependence 
on the mass of host clouds and then the effect of the ambient temperature. 

The values of swept gas mass, Mgw, at the fragmentation epoch are given in Table 1, 
which gives the minimum baryon masses of the host cloud necessary for the fragmentation 
under a given set of the ambient gas density (or collapse redshift) and the explosion energy. 
This mass varies by a factor of ~ 10 in the parameter range we computed from 6.07 x lO'' 
to 8.04 X 10^ Mq, and increases with the explosion energy and decreases for higher ambient 
gas density (or larger redshift at the collapse of host cloud). This means that less massive 
host clouds can form low-mass stars only when they are formed in the earlier epoch (higher 
z and thus larger po)- 

The necessary baryon mass, Mg^, of the host cloud for the fragmentation can be con- 
verted into the total (baryon plus dark matter) mass of the host clouds by assuming an 
appropriate model of the Universe. In Figure 9, we illustrate the permitted region (i.e., 
the baryon mass > Mg^) in the Eq — z plane for various total masses by assuming the 
Einstein-de-Sitter model, i.e., JIq = 1, f^s = 0.06, and Mp = Msyf{0,o/Qf,). For a given mass 
of the host cloud, there exists a lower limit to the ambient gas density or the redshift of the 
supernova explosion that allows fragmentation of the gas shell; the lower limit is an increase 
function of sq- In the case of the first collapsed objects with total masses of ~ 3 x 10^ Mq (see 
the hatched region of Fig.9), fragmentation can occur within the host cloud in the hatched 
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area of narrow parameters, e.g., if the supernova of Sq — 3x10^^ erg exploded at 2; 20 
or that of Sq = 10^"^ erg exploded at 2; ;^ 30. Because of the rather narrow range of Mg^, 
mentioned above, the mass of primordial clouds that can sustain supernova-triggered star 
formation is bounded sharply between the total mass of cloud Mp — 10® ~ 10^ Mq. That 
is, for the host clouds of Mt = 2 x lO^M©, supernova-triggered star formation is possible 
only in a narrow parameter region of 2; > 20 and £0 < 3 x 10^^ erg. For the host clouds of 
Mp > 1O^M0, star-formation can always be triggered, except for the narrow range in the 
right-bottom corner of z ~ 10 and £0 — 10^^ erg. 

We have also solved the thermal evolution of ambient gas consistently with the effect of 
the first generation stars taken into account, to show that the temperature in the ambient gas 
decreases efficiently to the range of The = 193 — 656 K by the time the fragment condition 
is satisfied (Table 1). By use of these temperatures, we may derive the condition that 
V > Cs,hc(3hc), which is plotted by a thick fine in the left bottom corner in Fig. 9. In the 
parameter region above and to the right of this line, i.e., if a supernova explodes with a larger 
explosion in the host cloud of higher density than delineated by this line, the gas shell swept 
by the SNR can fragment to give birth to low-mass stars surviving to date. Otherwise, the 
swept gas shell will dissolve and merge into the ambient gas without triggering fragmentation. 

In the above computations, we take into account the heating of ambient gas only by the 
SN progenitor, which ends with the supernova explosion. The temperatures in the ambient 
gas may vary and can be higher, however, if gas is irradiated by other heating sources, such 
as other first-generation massive stars or the SNR itself, and/or if the SNR shell expands 
beyond the Stromgren sphere, in which the temperature has not been raised high enough 
to ionize gas, and hence, the subsequent cooling is not so effective. In order to see the 
dependence of the evolution and fragmentation of the SNR gas shell on the thermal state 
of ambient gas, we calculate the evolution of several SNR models with the temperature of 
ambient gas set at a constant value, i.e.. The — ^o- For z — 20 and e — 10^^, the gas shell 
expands and cools enough to meet the fragment condition rg = Tdyn when the SNR radius 
reaches 92 and 85 pc in the models of constant ambient temperature The = 200 and 10^ K, 
respectively. The difference in the radii reflects the effect of deceleration by the ambient 
gas pressure. The condition under which the swept gas shell keeps expanding without being 
dissolved into the ambient medium {V > Cg^hc) is satisfied in the model of the lower ambient 
temperature, consistent with the above result of the calculation which solves the thermal 
evolution of ambient gas simultaneously. In the model of The = 10^ K, this condition is 
violated earlier, and hence, the gas shell has to be disturbed and merges into the ambient 
medium before it can gather gas sufficient to trigger the fragmentation. For such a high 
ambient temperature as The = 10^ K, the expansion of the SNR is halted at the radius of 
52 pc because of the high pressure of the ambient medium, and the swept gas shell will be 
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dissolved and melted in the host cloud. 

In Figure 9, we also show the parameter region, where the conditions of y = Cs,hc(^hc) are 
satisfied until the fragmentation epoch, by dotted lines for the different ambient temperatures 
of The = 200, 300, 500 and 10=^ K. 

In the parameter range left of these lines, where V < Cs(Thc), the SNR shell is dissolved 
into the ambient medium before the fragmentation epoch for a given ambient temperature 
The = ^0- If the ambient temperature is kept higher, the fragmentation demands a supernova 
of greater explosion energy, since the SNR has to sweep the amount of ambient gas necessary 
for fragmentation before the expansion is decelerated to sound velocity in the ambient gas. 
For example, if the ambient medium has been cooled to 300 K by radiative cooling, the 
explosion energy of 5o ^ 2 x lO^^erg is necessary for V > Cs,hc(300K) in a host cloud of 
Mp = 3 X 10^ Mq at z = 20, while for the cloud with The = lOOOK, a much stronger 
explosion energy oi Eq ^ 5 x lO^^erg is necessary for the SNR shell to fragment. For a high 
temperature of The ~ lO^K, such as assumed by Salvaterra, Ferrara, & Schneider (2003), 
SN-triggered star formation is possible only for a supernova (or hypernova) with an explosion 
energy exceeding 10^^ erg, consistent with their results. But this is not the case, since it is 
hkely that gas in the host cloud has been cooled after the supply of UV photons from the 
SN progenitor. 



5. Conclusions and Discussion 

In this paper, we have addressed two questions whether a primordial SNR can trigger 
star formation in a first collapsed object and whether low-mass stars can be formed to survive 
to date. For that purpose, we have studied the evolution of supernova remnants in primordial 
clouds with H2 and HD chemistry taken into account for a range of parameters of ambient 
density (or the redshift) of host clouds and explosion energy of supernovae, as summarized in 
Table 1. It is found that a gas shell swept by a supernova shock undergoes cooling and that 
HD molecules can effectively cool the temperature to as low as 32 — 154 K. Fragmentation 
begins when the mass in the gas shell reaches the range of 6.1 x 10'^ Mq — 8.0 x 10^ Mq. This 
is a necessary mass, Mgw, for star formation to be triggered by SNRs in primordial clouds of 
total mass Mt. Mg^ depends on the ambient density (or the redshift) and the SN explosion 
energy. The standard CDM model predicts that the masses of the collapsed objects from 
the 3 a peaks vary from Mp ^ 2 x 10^ Mq at z = 30 to ~ 10^ Mq at z = 10 (e.g., Haiman 
et al. 1996, Tegmark et al. 1997), which has a baryon mass of Mb ~ 10^ Mq {z = 30) and 
M5 ~ 5 X 10^ Mq (z = 10). If Ms„ < Mb = ^bMt, SNR can fragment within the host cloud 
to bring about supernova-triggered star formation. Otherwise, the gas shell passes the edge 
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of the host cloud, and hence, is blown out without fragmentation. 

In the range of parameters studied, fragmentation occurs in the pressure-driven expan- 
sion stage. The SNR expands to radii of 30 ~ 332 pc by this stage, which increases with 
the explosion energy and decreases with the ambient density. It takes 1.1 Myr — 45 Myr 
from the supernova explosion, much shorter than the Hubble age at that epoch, for the 
swept gas shell to start fragmentation, which is shorter for higher ambient density and for 
larger supernovae explosion energy. The lifetimes of massive stars of M ~ several tens Mq 
as r ~ lOMyr arc longer than the elapsed time for high z SNR {tf <^ 5Myr for z = 50 
models) and at least comparable for the supernova explosion at lower rcdshift (t/ — lOMyr 
for z = 20 models). Accordingly, it may be possible that fragmentation occurs even before 
other subsequent supernovae in the same primordial clouds contaminate the swept gas shell. 



We discuss the fate of host clouds after the low-mass stars have been formed in the gas 
shell. We evaluate the systematic expansion velocity V expected for newborn stars from 
the expanding velocity of a gas shell at fragmentation in Table 1. The escape velocities, 
V^sc) of the host clouds are also summarized in Table 1, and are evaluated from the total 
(baryon -|- dark-matter) mass contained in the SNR shell radius at the fragmentation epoch, 
Mt = Msw(JloM), as 



where Mgw represents the contained baryon mass. If V < V^sc? the second-generation low- 
mass stars remain within the host cloud. Since the mass of the host cloud is larger than 
the minimum mass of fragmentation Ms^{flo/flh), all the models calculated predict that 
newborn stars are bounded. 

Otherwise, the newborn stars will escape from the host cloud and be ejected into in- 
tergalactic space. This does not occur in our models. This tendency is consistent with the 
SPH simulation of SNR formed in low-mass minihalos M ~ IO^Mq by Bromm, Yoshida, 
& Hernquist (2003). That is, in a low energy model with Eq — lO^^erg the halo (cloud) is 
unchanged, while an SN with Eq — lO^^erg disrupts the halo (cloud). 

In our calculations, the fragments are likely to be bounded in the host cloud, because 
the velocity of the fragments is lower than the escape velocity estimated by the swept mass 
for all the models we studied. We only investigate the explosion energy of lO^^erg — 10^^ erg. 
If more energetic supernovae, such as pair instability supernovae with Eq > 10^^ erg occur. 



5.1. Fate of Stars Formed by This Mechanism 




(38) 
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newborn stars escape from the host cloud because their velocities are higher than the escape 
velocity. When the SNR escapes from the host cloud without fragmentation, the gas, which 
has large H2 and HD abundances, drifts in the intergalactic medium. 

5.2. Metal Abundance 

In the above, we have assumed the primordial abundances of elements for the gas shell. 
The fragments, and in particular, the stars formed out of them, are expected to be polluted 
by metals ejected from the SN, although the degree of metal abundance may vary with the 
efficiency of mixing between the SN cjccta and the surrounding material. Assuming complete 
mixing, we can evaluate the metal abundance in the gas shell. If the ejected mass of iron, 
Mej^Fe (Tsujimoto & Shigeyama 1998) is mixed with the swept mass, Mgw, this gives 

[Fe/H] ~ -2.8 + logio[(Mej,Fe/0.1 Me)/(M,w/5 x 10^ Mq)], (39) 

where the solar abundance of iron is assumed to be equal to 1.16 x 10~^. Figure 10 shows 
the expected value of [Fc/H]. From this figure, we found that the next generation stars 
triggered by primordial SNRs have a mctallicity of [Fe/H] ~ (~2.5) — (—4.5). The upper 
limit corresponds to low explosion energy sq = 10^^ erg and high mass iron ejection Mgj^Fe = 
0.5 Mq, while the other hmit corresponds to the high explosion energy £0 = 10^^ erg and 
low mass iron ejection Mej^pe = 0.05 M©. The actual distributions of metaUicity realized in 
the SNR, however, have to have a scatter around this mean value. It should be pointed 
out that the metal abundance of the low-mass stars formed by the supernova resembles that 
of the extremely metal-deficient stars in the Galactic halo [Fe/H]~ —3, discovered by the 
large-scale survey by Beers et al. (1992). However, the most metal deficient [Fe/H] = —5.3 
star that has been found, HE0107-5240 (Christlieb et al 2002), has too little metaUicity to be 
explained by the mixture of SN ejecta of a first-generation massive star and the primordial 
gas swept in its SNR shell. This seems to mean that HE0107-5240 was not formed in a cloud 
with Mr ~ several x 10*^ M0. 

How does the ejecta gas mix with a hundred times more massive cloud composed of 
primordial gas? Consider a cloud with smaller mass as Mp ~ 10^ Mq. The SNR shell 
sweeps all the gas in such a cloud, and the swept shell as well as the gas ejected from 
the SN are ejected into the intracloud space. In this case, metals mix not with the cloud 
medium, but with the intracloud medium. If this cloud belongs to a more massive density 
perturbation with Mt ~ 10^ Mq, such a system contracts later and forms stars with a 
hundred times more deficient metaUicity than the above-mentioned star formed in the cloud 
Mt several x 10^ Mq. 
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Specifically, our results suggest that the characteristic element abundances of ejecta from 
a single population III supernova can be imprinted on the abundance distributions in these 
extremely metal-poor stars. Their relevance is worthy of further study and if their connection 
is confirmed, we may gain information on the first collapsed objects in our Universe from 
the current status of these survivors. 

We have greatly benefited from discussion with A. Habe, R. Nishi, K. Wada, T. Mat- 
sumoto, K. Omukai, H. Uehara and T. Okamoto. Numerical calculations were carried 
out at the Astronomical Data Analysis Center, the National Astronomical Observatory of 
Japan. This work is supported in part by the Grants-in-Aid for Science Research from 
MEXT[09640308, 15204010 (MYF), 11640231, 14540233 (KT)]. 
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Table 1: Initial Conditions and Results 
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Symbols </, r, V, Vesc, Msv,, T, The, Cs,hc, [H2/HD] and [HD/H] represent the elapsed time, the radius of 
the shell, the expansion speed, the escape velocity, the swept mass, temperature of the shell, temperature 
of the ambient medium, sound speed of the ambient medium, and H2 and HD fractional abundances in the 
logarithmic scale, respectively. 
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A. Radiative Cooling Process 

As the radiative cooling process, we include the inverse Compton process owing to the 
cosmic background radiation, and line cooling of H, He, H2, and HD. Here in this Appendix, 
we briefly summarize the cooling rate. The unit of the cooling rate is erg cm^^ s^^ and that 
of the temperature is K. 



1. Inverse Compton cooling (Ikeuchi & Ostriker 1986) : 



Aie = 5.41 X 10-32 (1 + ^)^ i-f-] rie. 
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(Al) 



2. Helium cooling 

(a) CoUisional ionization cooling (Cen 1992) : 
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(b) Recombination cooling (Cen 1992) : 
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(c) CoUisional excitation coohng (Cen 1992) : 
AHe+,ex = 5.54 X 10-^^ T-'^-^gr 

AHe++,ex = 9.10 X 10-2^ 



/ rp \ 1/2' 
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3. Hydrogen cooling 



(a) CoUisional ionization cooling (Cen 1992) : 
Ah,c1 = 1.27 X 10~^^ 

(b) Recombination cooling (Cen 1992) : 
AH,e = 8.70 X 10-2^ (l^) " 

(c) CoUisional excitation cooling (Cen 1992): 

-1 

exp(-118348/r) uh- (All) 

4. Molecule Hydrogen cooling: Ahj is taken from the table of Flower et al. (2000). 

5. HD cooling: Ahd is taken from the table of Flower et al. (2000). 

6. Effect of the CMB radiation: 

AcMB = Ah2(2cmb) + Ahd(2cmb), (A12) 
where Tcmb is taken as 2.73(l+z). 
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B. Chemical Reactions 

We include the chemical reactions of 12 species: H, H+, H^, He, He+, He"'"^, H2, D, D+, 
HD, HD"*", and e~. In this Appendix we summarize the reactions we adopt. 



dt 



dnii+ 
dt 



driii- 
dt 

dt 
dniie+ 
dt 

(inHe++ 

dt 
dt 



dnji 
"df 

dnD+ 
dt 

driKD 
dt 
dniiD+ 
dt 



= k2 nH+ ne + 2 ki2 riHz rie + uh- ne + 2 ku riH- nH+ + k2o riB nH+ 

+ k22 nr, riHa + kg nn riR riR + 3 kio nm uh + ku ?T'H+ + 2 /cie JT-Hz '^Ha 

+ kiT riH - ki Ue — k-j nn rie - k^ nn JT-h- - ^21 ?T'D+ JT-h 

- k2A ?T'HD+ "-H - ^28 "-0+ "-H " 3A;9 riR riH JT-H - A;io riHa "-H 

- 2A;i5nHnHnH2 - 2A;i7nHnH, (Bl) 

= ki n-R He + /C21 + ^23 ^HD+ + /i^24 ^Ha + ^17 ^^H «H 

- /C2 nH+ rze - /Cii riH+ - /^u ?^H- ?^H+ - A;20 ?^H+ - k2e nuB wh+ 

- k27nT:,nji+, (B2) 

= /cr rie - kg riR rin- - ^13 n-u- - /cu riR- ?T'H+ , (B3) 

= ki nHe+ - fcs ^He ^e, (B4) 

= k^ riHe rie + nHe++ - ki nHe+ n-e - k^ nHe+ rie, (B5) 

^ k^ nHe+ rie - ^6 "'He++ "-e, (B6) 

= ksuu riu- + A;26 "^hd "h+ + ^i^g '^h "h ^^h + 2 /cis riu^ 

+ kie rin^ nn^ - hi nu+ - ku riu,, - k22 nn^ - ^24 ^0+ nn^ 

- kio nn - ki^ n^.^ - 2 kiQ 71^2 , (B7) 

= ki8 nD+ He + k2i nD+ riH + k25 riHD nn - kig no rie 

- k2o riB nH+ - k22 n-o nu^ (B8) 
= k2o riB nu+ + /c26 nuD nB+ - kis nr,+ rie - ^21 nD+ rin 

- k2A nD+ - k28 nD+ riu, (B9) 
= k22 riB riHz + ^23 ?T'HD+ J^H + ^^24 nji+ riHj - k2e riHD nii+ , (BIO) 

= /C27^D^H+ + fc28'T'D+ ^T'H - /E23«HD+ (BH) 



The reaction rates for the above reactions are taken from Palla, Salpeter, & Stahler 
(1983), Shapiro & Kang (1987), Abel et al. (1997), GaUi & Palla (1998) and Standi, Lepp, 
& Dalgarno (1998). 
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rate (cm ''s ^1 



reaction 



reference 



(1) H + e ^ H+ + 2 e h ^ exp[ 



(2) H+ + e ^ H + 7 



(3) He + e ^ He+ + 2 e 



32.71396786 
+ 13.536556xln(T/cv) 

- 5.73932875 xln(T/ev)2 
+ 1.56315498xln(T/ev)^ 

- 0.2877056xln(T/ev)^ 
+ 3.48255977 x 10-2xln(T/ev 

- 2.63197617 x 10-^xln(T/ev 
+ 1.11954395 X 10-^xln(T/cv 

- 2.03914985 x 10-6xln(T/ev) 

k2 = exp[ - 28.6130338 

- 0.72411256xln(T/cv) 

- 2.02604473 x 10-2xln(T/ev 

- 2.38086188 x lO-^xlnfT/ev 

- 3.21260521 X 10-^xln(T/ev 

- 1.42150291 X 10-5xln(T/ev 
+ 4.98910892 x lO-^xlnfT/ev 
+ 5.75561414 x 10-^xln(T/ev)^ 

- 1.85676704 x 10-«xln(T/evf 

- 3.07113524 x lO-^xlnfT/evf 

exp[ -44.09864886 

+ 23.91596563 xln(T/ev) 

- 10.7532302xln(T/ev)2 
+ 3.05803875xln(T/cv)3 

- 0.56851189xln(T/evf 
+ 6.79539123 x 10'^ x\n{T / evf 

- 5.00905610 X 10-''^xln(T/evf 
+ 2.06723616 x 10-^xln(T/ev)^ 

- 3.694916141 x 10-<^xln(T/ev) 



(4) He+ + e ^ He + 7 

(5) He+ + e ^ He++ + 2e 



k^ 
k5 



(6) He++ + e ^ He+ + 7 



(7) H + e ^ H- + 7 kr 



3.925 X 10 



-13 



(T/ev) 



-0.6353 



exp[ -68.71040990 

+ 43.93347633 xln(T/cv 

- 18.4806699xln(T/ev ' 
+ 4.70162649xln(T/ev)3 

- 0.76924663 xlnfr/evf 

+ 8.113042 X 10-2xln(T/ev)5 

- 5.32402063 x 10-3xln(T/ev)6 
+ 1.97570531 x IQ-^xlnfT/ev)^ 

- 3.16558106 x IQ-^xlnfr/evf 



k, = 3.3610-0 T-V^ (^) 



-0.2 



0.7 



1.0 X 10-18 T 



dex[ - 14.10 

+ 0.1175 log T 

- 9.813 X 10-=^ (logr)2 ] 



T < 1.5 X 10^ K 
T > 1.5 X 10^ K 
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reaction 



rate ( 



cm 



reference 



(8) H + H- ^ H2 + e 



(9) H + H + H ^ H2 + H 

(10) H2 + H^H + H + H 

(11) Ho +H+ ^ H+ + H 



(12) 
(13) 



H2 + e ^ 2H + e 
H- + e ^ H + 2e 



(14) 
(15) 
(16) 

(17) 
(18) 
(19) 
(20) 



exp[ - 20.06913897 T > 1160 K 
+ 0.22898xln(T/ev) 
+ 3.5998377 X 10-2xln(T/ev)2 

- 4.55512 X 10-3xln(T/ev)'^ 

- 3.10511544 X 10-^xln(T/ev)^ 
+ 1.0732940 X 10-^xln(T/ev)^ 

- 8.36671960 x 10'^ xln{T/ evf 
+ 2.23830623 x lO-^xlnfr/ev/ 

1.428 X 10-9 



T < 1160 K 



kg = 5.5 X 10 

ho 
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X T- 



k 



11 



ki2 
ki3 



= 6.5 X lO-^T-1/2 exp(-52000/r) 
x[l-exp(-6000/T)] 

= exp[- 24.24914687 

+ 3.4008244 x ln(T/ev) 

- 3.89800396 x ln(T/cv)2 
+ 2.04558782 x ln(T/ev)3 

- 0.541618285 x ln(T/ev)^ 

+ 8.41077503 x 10"^ x ln(r/ev)^ 

- 7.87902615 x 10"^ x ln(T/ev)6 
+ 4.13839842 x 10"^ x ln(T/evj^ 

- 9.3634588 x 10"^ x ln(T/ev)^ ] 

= 5.6 x 10-^^rV2 exp(-102124/T) 

=exp[- 18.01849334 

+ 2.3608522 x ln(T/ev) 

- 0.28274430 x ln(T/ev)2 

+ 1.62331664 x lO^^ x ln{T/evf 

- 3.36501203 x lO'^ 
+ 1.17832978 x lO^^ 

- 1.65619470 X lO^^ 
+ 1.06827520 x 10"^ 

- 2.63128581 x 10"^ 



X ln(T/ev)4 
X ln(T/evf 
X ln{T/ev f 
X ln(r/ev)^ 
X ln(T/evf 



H- + H+ ^ 2H 

H + H + H2 ^ H2 + H2 

H2 + H2 ^ H + H + H2 

H + H^H+ + e + H 
D+ + e ^ D + 7 
D + e ^ D" + 7 
D + H+ ^ D+ + H 



k 



-1/2 



14 



(100) 

-30rp-l 



= 7 X 10- 
ki5 = 6.875 X 10- 

ki6 = 8.0 X 10-8^-1/2 exp(-52000/T) 
X [1 - exp(-6000/r)] 

kn = 9.86 X 10-i^r°-5exp(-158000/T) 
kis = 3.6 X 10-12(2^/300) 
kig = 3.0 X 10-i6(r/300)°-95exp(-r/9320) 
k2o = 3.7 X 10-iorO-28gxp(_43/j.) 



3 
3 

1 
1 



1 
3 
3 

3 
4 
4 
4 
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reaction rate (cm '^s ^) reference 

(21) D+ + H ^ D + H+ k2i^ 3.7 x 10-^° T^ ^s 4 

(22) D + H2 ^ H + HD A;22 = 9 X IQ-^^ exp(-3876/r) 4 

(23) HD+ + H ^ H+ + HD A;23 = 6.4 x 10-^° 4 

(24) D+ + H2 ^ H+ + HD k24 = 2.1 x lO"'^ 4 

(25) HD + H ^ H2 + D A;25 = 3.2 x IQ-^^ exp(-3624/r) 4 

(26) HD + H+ ^ H2 + D+ k2e = 1.0 x IQ-^ exp(-464/r) 4 

(27) D + H+ ^ HD+ + 7 k27^ dex [- 19.38 - 1.523 x logT 5 

+ 1.118 x (logT)2 
- 0.1269 X (logT)3 ] 

(28) D+ + H ^ HD+ + 7 k28 =dex[- 19.38 - 1.523 x logT 5 
+ 1.118 X (logTy - 0.1269 X (logTf] 

Table 2: References.— 1. Abel et al. (1997); 2. Shapiro & Kang (1987); 3. Palla, Salpeter, 

& Stabler (1983); 4. GaUi & Palla (1998); 5. Standi, Lepp, & Dalgarno (1998). 
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Fig. 1. — Evolution of SNR radius (top panels) and the shell expansion velocity (bottom 
panels) is compared among the models of different ambient densities {z=10, 20, 30, 40, and 
50) but a fixed explosion energy {eq = 10^^ erg; left panels) and among the models of different 
explosion energies {eq = 1,3,5, andlO x lO^^erg) but a fixed ambient density {z = 20; right 
panels). 




Fig. 2. — Time variations in the temperature of the gas shell (the top and middle panels) 
and that of the ambient gas (the bottom panel) are plotted against the elapsed time. In 
the top and bottom panels, models with different ambient densities (or different redshifts of 
formation of host cloud z = 10, 20, 30, 40, and 50), but a fixed explosion energy {eq = 10^^ 
erg), are compared. The middle panel compares the models with different explosion energies 
{eq = 1, 3, 5, and 10 x 10^^ erg), but a fixed ambient density (or redshift z = 20). 
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Fig. 3. — Variations of the pressures of the swept gas, P, the cavity gas, Pin, the post-shock 
gas Pout, and the ambient gas, Phc- Model A20 of {z,eo) = (20, 10^^ erg) are shown for the 
ante-fragmentation phase. 
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Fig. 4. — Time variations in the number density of gas in the gas shell for models with 
different energies {eq = 1, 3, 5, and 10 x 10^^ erg), but a fixed redshift of z = 20. 
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Fig. 5. — Evolutions of four typical timescales {rg: free-fall timescale, r^jn'- dynamical 
timescale, Texp: expansion timescale, and Tcooi: cooling time scale) are plotted against the 
elapsed time. The left panel is for a model with {z, eo) = {20, Ix 10^^ erg), and top-right panel 
is for a model with (z, eo) = (50, 1x10^^ erg) and the bottom-right panel is for a model with 
{z,eo) = {20, 5x10^^ erg). Texp is shown only for the Sedov- Taylor stage. 
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Fig. 6. — Time variations in the fractional abundances for 10 species (H, H"*", H~, H2, He, 
He^, He^+, HD, D, and D"^) in the gas shell for the model with Eq = 10^^ erg and z = 20. 
The abundance is plotted as a function of temperature (left panel) before the fragmentation, 
and as a function of number density (right panel) after the fragmentation. In the post- 
fragmentation phase, we assume a cylindrical fragment contracts in the radial direction. 
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Fig. 7. — Cooling rates accomplished in the gas shell (Aaii: total cooling rate, Ac: inverse 
Compton cooling, Ah: Hydrogen coohng, Aue. Helium cooling, A^^: molecular Hydrogen 
cooling, and Ahd^ HD molecular cooling) are plotted against the gas temperature in the ante- 
fragmentation phase (left panel) and against the number density in the post-fragmentation 
phase (right two panels) . The top-right panel indicates the cooling rates in the case of spher- 
ical collapse, and the bottom right is for the case of cylindrical collapse. This corresponds 
to the model with {z,eo) = (20, lO^^erg). 
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Fig. 8. — Evolution of gas temperature (T: solid curves) and the Jeans mass (Mj: broken 
curves) in the fragments during the post-fragmentation phase for a model with (z, Eq) = 
(20,10^^ erg). Two geometrical models of collapse are shown for spherical (Tg, Mjs) and 
cylindrical (Tc, Mjc) collapses. 
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Fig. 9. — Lower bounds of the redshift z (the left axis) or of the background density, po^ in the 
host clouds (the right axis) for star formation to be triggered by a single supernova explosion 
are plotted as a function of the explosion energy, of supernova for the host clouds of total 
mass, given mass Mt= 2, 3, 5, and 10 x 10^ Mq, from the top to the bottom. The thick line 
denotes the border oiV = Cs,hc(7'), which divides whether the gas shell swept by SNR can 
fragment or merge into the ambient gas, with the temperature in the ambient gas derived 
by solving the chemical reaction and thermal evolution of the ambient gas simultaneously. 
E.g., the shaded region denotes the parameter range in which we expect the SNR-induced 
fragmentation of low-mass second-generation stars in the host clouds with the total mass of 
Mr = 3 X IO^Mq. Dotted lines represent the dependence on the shell expansion speed (or 
the temperature)in the ambient gas for the temperatures of 200 K, 300 K, 500 K, and 1000 
K. They give the boundaries for the gas shell, swept by SNR, to fragment without mixing 
into the ambient matter, i.e., supernova-triggered, low-mass star formation is expected in 
the parameter range to the right of these lines, and requires a supernova of larger explosion 
energy for higher temperatures in the ambient gas. For the supernova of weaker energy, the 
SNR shell dissolves into the ambient gas before the fragmentation. 
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Fig. 10. — The expected metallicity of [Fe/H] is plotted against the initial z. [Fe/H] is 
estimated by the ratio of the swept mass and the ejected iron mass. Ejected iron mass is 
assumed to be MFe,ej = 0.1 (thick center lines), 0.05 Mq (left border lines) and 0.5 Mq 
(right border lines) estimated by Tsujimoto & Shigeyama (1998). Solid and dotted lines 
represent the model of eo = 10^^ erg and that of eo = 10^^ erg, respectively. 



